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Abstract 

In this article we study the positive solutions of the parabolic semilinear system of 
competitive type 

j u t - Au + v p = 0, 
\ v t - Av + u q = 0, 

in il x (0, T), where fl is a domain of R N , and p, q > 0, pq ^ 1. Despite of the lack of 
comparison principles, we prove local upper estimates in the superlinear case pq > 1 of 
the form 

u(x, t) <: Ct- (p+1)/(p ' ? - 1) , v{x, t) ^ C<- (9+1)/(f " ? - 1) 

in u) x (0,Ti) , for any domain wcc!! and T x € (0,T) , and C = C(N,p,q t Ti,u). For 
p, q > 1, we prove the existence of an initial trace at time 0, which is a Borel measure 
on n. Finally we prove that the punctual singularities at time are removable when 
p, q ^ 1 + 2/N. 
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1 Introduction 



Let Q be a domain of 1^ (N > 1) and < T ^ oo. In this work we are concerned with the 
positive solutions of the parabolic system with absorption terms 

ut - Au + V p = 0, 

v t - Av + u q = 0, ^ ' ' 

in f2 x (0, T) , with p,q > 0, pq ^ 1, in particular in the superlinear case where pq > 1. 

This system appears as a simple model of competition between two species, where the 
increase of the population of one of them reduces the growth rate of the other. Independently 
of the biological applications, it presents an evident interest, since it is the direct extension 
of the scalar equation 

U t - AU + U Q = 0, (1.2) 

with Q ^ 1. For Q > 1, any nonnegative subsolution of equation (jl.2p in Q x (0, T) satisfies 
the following upper estimate: for any bounded C 2 domain uj C O 

U(x,t) S({Q -l)tr 1/{Q ' 1] +Cd{x,duj)- 2/ ^ Q - 1) V(i,t)Gwx(0,r), (1.3) 

where d(x, du) is the distance from x to the boundary of u and C = C(N, Q), see |15| . This 
estimate follows from the comparison principle, as shown at Proposition 13.41 Moreover it 
was proved in |15| that any solution U of equation (jl.2p in (7 x (0, T) admits a trace at time 
in the following sense: 

There exist two disjoints sets 1Z and S such that 1ZU S = Q, and TZ is open, and a 
nonnegative Radon measure [x on TZ, such that 

• For any xo € TZ, and any ip G C®{TZ), 

lim / U(.,t)iP= f $dn, 



For any open set U such that hi R 5 ^ 



lim / u(.,t) = oo. 

Moreover the trace (S,fj,) is unique whenever Q < 1 + 2/iV. 

Up to now, system (jl.ip has been barely touched on. Indeed an essential difficulty 
appears: the lack of results for comparison principles. As a consequence, most of the 
classical properties of equation (II. 2j) . based on the use of standard supersolutions, are 
hardly extendable. Some existence results are given in [12] for bounded initial data, and 
then in [3] for more general multipower systems with non smooth data, see also [13] for 
quasilinear operators. Otherwise the existence of traveling waves is treated in [9]. For the 
associated elliptic system 

-Au + t* = 0, , , 

-Av + u« = 0, 1 ' ; 
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the isolated singularities are completely described in [3] for the superlinear case pq > 1 
and for the sublinear case pq < 1, see also [T7], [18] for p, q ^ 1. The study shows a great 
complexity of the possible singularities; in particular many nonradial singular solutions 
are constructed by bifurcation methods. The question of large solutions of system (jl.4p is 
studied in the radial case in [TT], showing unexpected multiplicity results, and the behavior 
of the solutions near the boundary is open in dimension N > 1; the existence is an open 
problem in the general case. For such competitive problems, some more adapted sub- 
supersolutions and super-subsolutions have been introduced, see [H], [16], [3], [10], but the 
problem remains to construct them. The uniqueness is also a difficult problem, as it was 
first observed in pQ. 

Our first result consists in local backward upper estimates for the solutions of the system: 
defining the two exponents 

a= -, b = -, (1.5) 

pq — 1 pq — 1 

we obtain the following: 

Theorem 1.1 Assume that pq > 1. Let (u,v) be a positive solution of system M-l}) in 
Q x (0, T) . Then for any domain wCCSi (o> = R N if Tl = R N ), 

u(x, t) £ Ct- a , v(x, t) £ Ct~ b , V(x, t) e uj x (0, T) , (1.6) 

for some C = C(N,p,q,T,u)). 

Our second result is the existence of a trace in the following sense: 

Theorem 1.2 Assume that p,q > 1. Let (u,v) be a positive solution of the system in 
x (0, T) . Then there exist two disjoints sets 1Z and S such that 7Z U S = f2 ; and 1Z is 
open, and nonnegative Radon measures [ii,[i2 on T^> such that the following holds: 

• For any xq € 1Z, and any t/> € C®(1Z), 

lim / u(.,t)i/j = / tpdfii, lim / v(.,t)ip= / tpdfj,2- (1-7) 
t^Jn J n t-+°Jn Jn 

• For any open set U such that U PI S 7^ 0, 

lim / (u(.,t) + v(.,t)) = 00. (l.i 



As a consequence we can define a notion of trace of (it, v) at time 0: 
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Definition 1.3 The couple 8 = (81,82) of Borel measures 81,82 on £2 associated to the 
triplet (S, m, H2) defined for i = 1,2 by 

m(E) if ECU, 

l[ ' ~ \ 00 if Ens 7^0, 

is called the initial trace of (u,v). 

Finally we give a result of removability of the initial singularities inspired by [6] Theorem 

2]: 

Theorem 1.4 Assume thatp,q ^ 1+2/N. If there exists a positive solution (u,v) of system 
(COP infix (0, T) such that 

lim/ (u(.,t)+v(t))<p = 0, V^C c °°(n\{0}), (1.9) 

then u, v G C 2,1 (0 x [0, T)) and u(x, 0) = v(x, 0) = 0, Vx G O. 

In each section we point out some questions which remain open. 

2 Some existence results 

Next we recall some results that we obtained in [3] where we studied the existence and the 
eventual uniqueness of signed solutions of the Cauchy problem with initial data (uo,^o) 

u t — An + \v\ p H _1 u = 0, . . 

A 1 I 19 I 1-1 n 

vt — A« + \u\ \v\ v = 0, 
where p, g > 0, and 

( 1 if u > 0, 

Inp 1 n = < if u = 0, 
( -1 if u< 0. 

In particular we showed in [3] the following results: 

Theorem 2.1 Assume that Q is bounded. Suppose that uq G L e (Q) and vq G L x (U) with 
1 ^ #, A ^ 00, with 

max(^) < 1 + 2/N, 

or that uo,vo are two bounded Radon measures in ft, and 

m&x{p,q) < 1 + 2/N. (2.2) 

Then there exists a weak solution (u, v) of the system with Dirichlet or Neuman conditions 
on the lateral boundary, such that for any ip G C^(Q), 

lim / u(.,t)ib = / i^duo, lim / v(.,t)ip = / ij)dvo. (2-3) 
t->oJn Jn t^Jn Jn 
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Also, there exist two solutions (u\,vi) and (^2,1*2) such that any solution (u,v) satisfies 
U\ ^ u ^ U2 and V2 S v = v i- 

Moreover, ifp,q > 1 and uo G L e {Q) and v$ G L A (f2) with 

i/ien (u, -u) is unique; in particular this holds for any uq,vq G L^fi), i/ h2. '3) is satisfied, or 
ifu ,v G L e (ft) wi/j > N(m&x(p,q) - l)/2. 

3 Local a priori estimates 

When looking for local upper estimates of the nonnegative solutions of system (|2.ip near 
t = 0, we notice that the system admits the solution (0, v ) with v a solution of the heat 
equation in Q x (0,T), for which we have no estimate, since the set of solutions is a vector 
space. That is why we suppose that u and v are positive in £1 x (0, T) . The question of upper 
estimates for one of the functions is very closely linked to the question of lower estimates 
for the other one. 

We define a solution of problem (jl.ip in fix (0, T) as a couple (u, v) of positive functions 
such that u G £ 9 OC (0 x (0,T)), u G £f oc (^ x (0,T)) and 

(-utpt - uAip + v p cp) = 0, (3.1) 

f2x(o,r) 

(-wpt - wAp + «V) = 0, (3.2) 

Qx(0,T) 

for any 99 G P(0 x (0, T)). From the standard regularity theory for the heat equation it 
follows that u,v G C/ ' c (^ x (0, T)), and then u, v G C°°(0 x (0,T)) since u, u are positive. 

As in the case of the scalar equation (jl.2p , the system (jl.ip admits a particular solution 
(u*,t>*) for pq > 1, defined by 

u*(i) = A*t~ a , v*{t)=B*t~\ 

where 

= (p + + i)P(pg _ l)"(P+l), (B*) pq ^ 1 = (q+ l)(p + l)' (pg - l)-( 9+1 ). 



In [I], the authors studied the singularities near of the positive solutions of the asso- 
ciated elliptic system (|1.4p in B(0, 1)\{0}. System (jl.4p admits particular solutions when 
min(2a, 2b) > N — 2, given by 

w*(x) = ^4* |x| , f*(x) = -Z?* |x| , 

with 

API- 1 = 2a(2a + 2- N)((2b(2b + 2- N)f, B™- 1 = 26(26 + 2 - iV)((2o(2a + 2 - N)f. 
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When pq > 1 the following upper estimates hold near near : 

u{x) < C \x\~ 2a , v(x) < C |xp 26 , 

for some C = C(p, q, N). The proofs were based on estimates of the mean value of u and 
v on the sphere {|x| = r}, on the mean value inequality for subharmonic functions, and a 
bootstrap technique for comparisons between different spheres. 

For system (jl.ip the estimates (|1.6p are based on local integral estimates of the solu- 
tions, following some ideas of [5] for elliptic systems with source terms. Then we use two 
arguments: the mean value inequality in suitable cylinders for subsolutions of the heat 
equation, and an adaptation of the bootstrap technique of 



Notation 3.1 For any cylinder Q = u x (s, t) C 0, x (0,T) and any w G L l (Q) we set 



1 

w = 



1 [ 
/ w. 

S J bj 



Q 

For any p > 0, we define the open ball B p = B(0,p) and the cylinder 

Q p = B p x [-p 2 ,0] . 

We denote by £i the first eigenfunction of the Laplacian in B\, such that J B £i = 1, with 
eigenvalue Ai, and by £ the first eigenfunction in B p with eigenvalue Ai jP = Ai//> 2 , defined 
by 

£(*)=&(-), VxeSp. (3.3) 
First we need a precise version of the mean value inequality. 

Lemma 3.2 Lei be any domain in M. N , and let w be a subsolution of the heat equation 
in Q x (0, T), with w E C 2,1 (Q x (0,T)). Then for any r > 0, there exists a constant 
C = C(N,r), such that for any (xo,io) an< ^ /° > such that (xo,to) + Q P Cflx (0,T), and 
for any e G (0, 1/2), 

sup u> ^ Ce'^ (-1 _ w r J . (3.4) 

Proof. This Lemma is given in case e = 1 in [8] for solutions of the heat equation, and we 
adapt its proof with the parameter e. We can assume that (xq, to) = and r € (0, 1) .From 
[8] there exists Cjy = C(N) > such that for any a G (0, 1) , 

sup w S C N (1 -cr)" (iV+2) | w. (3.5) 



G 



For any n € N, let p n = p{\ - e)(l + e/2 + ... + (e/2) n ), and M n = supg \w\ . From ([33 
we obtain 

Mn s c N (i - -Ed-)-(n+h I w . 

Pn+l J Q Pn+1 

thus with a new constant Cn 

M n ^ C^- (n+1)(7V+2) / w. 
From Young inequality, for any 6 € (0, 1), we obtain 



M n <,C N e-^ N + 2 )M^ x jf_ 



T 

w 



^ W n+1 +r5 1 -^(C iY ^ (n+1)(W+2) )^ 
Defining D = rJ 1 - 1 ^^ and 6 = £ -(^+ 2 )A, we find 

M„ £ 5M n+ i + b n+1 D ^ £ 
Taking 5 = 1/26 and iterating, we obtain 

M = sup \w\ S S n+l M n+ i + bDJ2(6by I I 
S d n+1 M n+l + 2bD 



w 



QPn + 1 



w r 



p(l-e) 1 = \ 

1 



W 

'Pn + l / 

Since Q Pn+1 C Qp(i+e)i we deduce (|3.4p by going to the limit as n — > oo. ■ 
Next we recall a bootstrap result given from [2J Lemma 2.2]: 

Lemma 3.3 Let d,h,£ S R iirai/i d € (0,1) and y, be two continuous positive functions 
on some interval (0, R] . Assume that there exist some C ,M > and Eq £ (0, 1/2] such that, 
for any e G (0,£ ], 

y(r) S C e~ h $(r) y d [r(l - e)] and max <&(t) ^ M $(r), 

re[r/2,r] 

or eke 

y(r) ^ C e _/l $(r) y d [r(l + e)] and max $(r) 1^ M $>(r), 

re[r,3r/2] 

/or any r £ (0, i?/2] . T/ien i/iere exists another C > suc/i i/iai 

y(r) ^ C ^(r) 1 ^ 1 -^ 

on (0,R/2] . 
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Next we prove the estimates (|1.6p . 

Proof of Theorem 11.11 We consider any point (xo,to) G f2 x (0, T) , and any p > 
such that B(xq,p) = xq + B p C By translation we can reduce to the case xq = 0. For 
given s S (0,1), we consider a smooth function r]o(t) on [— 2s,0] with values in [0,1] such 
that r]Q = 1 in [— s, 0] and r/o(— 2s) = and (vo)t(t) 5j Cs _1 . Choosing s such that 
< to — 2s < to, we set r](t) = rjo(t — to). We multiply the first equation in (jl.lj) by 

where £ is defined at (|3.3|) . and A > 1, which will be chosen large enough. We obtain 

j\! <VWJ + / ^V = A/ <V~V(*)+ / «(A^ A )r/ A . (3.6) 



By computation, we find 



P 2 A£ A (*) = ac a (-; 



-AA 1 e A (-) + A(A-lK A - 2 |V6| 2 (-)- 

p p 

x/e'-2 



For given ^ > 1, if A > 2£', the function gt(y) = is bounded, thus 



u (.,t)(^V(^ 



B„ 



< 



A(A - 


1) 


P 2 




A(A - 


1) 


P 2 




A(A - 


1) 



u(x,t)(^- 2 \V^\ 2 )(-) V x (t)dx 
u{x,t)^gi{-)r, X {t)dx 



l/t 



and even with different constants C = C (N, £) 



[ U (.,t)\Ae\v x (t)sx\ lP - 2 f U (.,t)ev x (t)+cp N / £ '- 2 ( [ 

J Bp J Bp \ J Bp 



U (.,tyev x (t) 



III 



l/i 



(3.7) 



Moreover 



1/f 



B 



SCp N ^' s -U / M (.,i)WW 



s 



Integrating (|3.6j) on (to — 2s, to) > and using Holder inequality, 



JB D 



+ 



to 



t -2s JB : 



?er} X ^Cp N l l \p- 2 + s- 1 ) 



t -2s \JB. 



l/i 



^Cp N ^{p- 2 + s- l )s 1 ^ 



In the same way, for any k > 1, if A > 2k', 

l-to 



[ v(.,t )e 

Jb 



^W^C P N / K '(p- 2 + s- 1 )s l ^'( 



I to -2s J B p 

Next we discuss according to the values of p and q. 



to 



to -2s J Bp 



to 



t -2s 



l/i 



u i e-n x I • (3-8) 



• k £V) 1/k - (3-9) 



First case: p,q > 1. We take £ = q,K = p, and 2s = /? 2 and consider any to such that 
< to — p 2 < to < T. Let us denote Q p = (0, to) + Q p - Then 



/ 1 S Cp( N+2 ^'- 2 (if 

J J Q n \JJQ, 



U q £ X T) X 



1/9 



l/p 



that means 



/ ^ A r/ A S Cp- 2 ( f u q i 
J) Q P \J Q P 

I u q ev x scp- 2 [J v n 

J) Q P \J Qp 



v 



1/9, 



1/p 



(3.10) 



Hence 



Then we get an estimate of the form 

u q ) 1/q S 



l/pq 



V 



c 



and similarly 



Qp/2 

( / v p ) l / p < 



p 2(p+l)/(pq-l) 
C 



(3.11) 



q p/2 fh+Wte-i) 

But u is a subsolution of the heat equation, hence there exists a C = C(N,q) such that 

u(x,t)SC(Jf u q )V q , 

J> Qp/2 



(3.12) 
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from Lemma 13.21 with r = q and e = 1. Taking p 2 = to/2M, with M > 1, we deduce the 
estimates 

U{ - X > > = t (p+l)/(p</-l) ' = t (9+l)/(p5-l) ' 

for any t G (0,T) and any x G such that £(x, y / t/2M) C O, with C = C(JV,p, g, M).Then 
L6l) follows. 



General case: > 1. We can assume p ^ 1 < q. Taking again < to — p 2 < to < T and 
2 s = p 2 , and using (|3.8p with I = q > 1, we find again (|3. lOj) . Using (|3.9p . we find for any 

k > 1, 

(\ 1 / K / \ 1 / K 

/ [ ^ V s c^ n+2 ^ k '- 2 sup v 1 -?/* 
J JQ P J Q P 

(3.13) 

More precisely, for any e G (0, 1/2), from Lemma 13.21 we find taking r = p and k = q, 



(\ Vp 
// M 
^Qp(l+e) / 



then 



Qp \JJQp J \jJQp(l+e) J 

/ \ Vp 

= C e-^^p-^ 2 ^ ff A 



\J JQp(l +e ) / 

Using (|3.13p we deduce 

u« S Ce-^^^p^ 2 ^'- 2 -^ 2 ^ iff vA ^ 



p(l-e) V ^P(l+e) / 

= fo-t^+l^^l^WWl/ll-l/p)^ iff v p\ . 

\JjQp(l+e) ) 

setting h = 2X + (N + 2)(q — p)/p 2 q, that means 

(\ Vp 
/ v p \ . (3.14) 
Qp(l+e) / 

Next from (|3.10p we have 

ff v^Cp- 2 (ff J] , (3.15) 
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thus changing p(l — e) into p(l + e), 

1/9 



■» Qp(l+e) y" Qp(l+4e) , 



Hence from (|3.14p . we deduce 

u q S Ce- h p- 2 ^/ p (f u«Y . 
Qp(i-e) y" Qp(i+4 £ ) / 

From Lemma 13.31 we conclude that 

\ (P9-l)/9 

Hence (|3.11|) follows as above, and then (|3.12|) from (|3.15|) . and the conclusion follows again. 
■ 

Next we give a first extension of the scalar estimate (|1.3|) to system (jl.ip . using some 
ideas of 0J p. 243]. 

Proposition 3.4 Let q ^ p > 1. Let (u,v) be any positive solution of system (COP * n 
17 x (0, T) , where ft is a bounded C 2 domain Then there exists a constant C = C(N,p,q) 
such that 

u^ +i y ( P +1 \x,t) + v(x,t) ^ C{t + d 2 {x,dn))- 1/(p - 1 \ V(x,t)£flx (0,T) (3.16) 

Proof. Let F = (k + u) d + u, with d = (<? + l)/(p + 1) > 1 and fc > 0. Then 

F t - AF = d(k + u) d - l (u t - An) - d(d - l)(fc + u) d ~ 2 [Vu| 2 + «t - Av 
S -d{k + u) d ~ 1 v p -u q . 

But (k + u) q ^ 2i- 1 (ki + u q ), thus 

F t -AF + d(k + u) d ~V + 2 1 " 9 (A: + u) q ^ k q . 

Observe that (k + u) q = (k + u)^ 1 ^ + it)*, and F p ^ 2P" 1 ((/c + «)* + «p). Then 

F t - AF + c(k + u) d - l F p S k q , 

with c = 2 1_p min(d, 2 1_9 ). In particular, taking = c~ 1 ^ d ^ 1 \ F is a subsolution of equation 

U t -AU + U p = K (3.17) 

in x (0,T) , where K = k q = K(p,q). Let /(i) = ((p - l)t))~ l/{Q ~ X) and let 5 be the 
maximal solution of the stationary problem —AU + U p = in fi such that g = oo on <9f2. 
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Then for any e > 0, the function (x,t) h4 G £ (x,t) = K x l p + / (t — e) + g{x) is a supersolution 
of equation ()3.17j) in fi x (e, T) . Going to the limit as e — > 0, it follows that 

F{x,t)ZKV* + f(t) + g(x) 

infix (0, T) ; then there exists a constants C = C'(N,p) such that 

F(x, t) S KVp + f(t) + C'd{x, 3fi)- 2 /(p-i) , V(x, t) € fi x (0, T) , 

and the conclusion follows. ■ 

Open problem: The estimate (|3.16p does not appear to be optimal, except in the 
case p = q where u = v is a solution of the scalar equation (|1.2|) . Can we obtain for p, q > 1, 
and even for pq > 1, an upper estimate in fi x (0, T) of the form 

u(x, t) S C(t + d 2 (x, dn))~ a , v(x, t) ^ C(t + d 2 (x, dfl)y b , 

with C = C{N,p,q)l 



4 Initial trace 

First we show some properties available for any p, q > 0. 



Lemma 4.1 Assume p, q > 0. Let (u,v) be any positive solution of system and let 

B(x , p) C fi. If Jq f B (x ,p) yP < °°> then $B{x ,p) u i-^) is bounded as t -> for any p < p, 
and there exists a Radon measure mi tP on B(xq,p) such that for any ip G C£°(B(xq, p)), 

lim / u(.,t)%[) = mi, p (ip). 



>B(x ,p) 

Proof. We reduce to the case xq = 0. We set 



X(t)= u(.,t)t, Y(t)= v(.,t)t, Z(t)= ui(.,t)e, W(t)= vf(.,t)t. 

<j Bp J Bp J Bp J Bp 

(4.1) 

where £ is defined at (|3.3p and A ^ 2. We obtain 

X t + W=j([ u( x ) + [ v^ x = / u(A^ 
= -AAi 



Bp ) J B a J B D 



hence 



, p [ < A + A(A-1) ! < A ~ 2 |V£| 5 

J Bp J Bp 

^ -AAi, p f < A = -AA liP X, 

JB 
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By integration we obtain for any t < 9 

e XXl ^ e X{9) - e XXl ^X(t) + J e XXl ^ s W{s)ds ^ 0; 

and from our assumption W € L 1 ((0,T)). Then e XXl ' pt X(t) is bounded, and in turn 
X(t) is bounded. Then f B u(.,t)^ x is bounded, hence f B r Xo p )U(.,t) is bounded. Let 
ip e C™(B(xo,p)). Then 

J Bp J Bp */ Bp 

Since Aip is bounded with compact support, we have + \Atp\ ^ C£ A for some positive 
constant C, and thus J B u(Aip) is bounded, implying j B u(., t)ip has a limit mi :P (ip), which 
defines a Radon measure rni )P on B p . ■ 

Lemma 4.2 Assume p,q > 0. Xei (n, u) 6e any positive solution of system hl.l}) , and let 
B(xo,po) C Q. If f B r X0 po \ u(.,t) is bounded as t 0, t/ien 

(%) for any p < po, J. J B v p is bounded; 

(ii) for any p < po, any 1 f$ a < 1 + 2/iV, and any < t < 9 < T 

[ [ u a dx ^ C, (4.2) 

where C = C(N,p,q,p,po,a). 



Proof. We still reduce to the case xq = 0. 

(i) Let < t < 9 < T with fixed 9, and C = sup {ofi] j B u(.,t). Let ip G C™(B po ) with 
values in [0, 1] such that t/j = 1 on B p . Taking ip as a test function in the equation in u and 
integrating between t and 0, we find 

u(Atp) S C\\Aip\\ Loo( Q^ , 
thus 




u(.,9)i>+ / ^ C(ll AV|lL» ( n) + 1), 

Jt Jb po 

hence f t J B v p is bounded. 

(ii) Here we use the ideas of [2, Propositions 2.1,2.2.] relative to quasilinear equations 
in order to estimate the gradient. Since a < 1 + 2/N, we can fix a = a(o~) such that 

- 1< a < and a S " + 1 + 2/N. (4.3) 
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Let p be fixed such that p < p < po. We multiply the equation in u by (1 + u) a ^ x , where £ 
is defined at (|3.3p . with A ^ 2/ |a| . Then we find for fixed 9 < T, and any < t ^ 9 



l — I (i +u (.,t)r +i e+\<x\ I I (l+ur-^vufe 

— J-r / (l + u(.,9)) a+l C X + f I v p {l + u) a i x + \f f {l + u) a i x - l Vu.Vi. 



Applying twice the Holder inequality, we find 

~~r I {l + u{.,t)T +i e + \\a\ II {l + uY~ l \Vu\ 2 e 

+ 1 J B 1 Jt JB 

v p (l + u) a £ x + C 



a + 
<C + 




t JBn 




(i+ u ) Q+i £ A - 2 |ve 



t J B„ 



<c + 



\a\ 



X-2/\a\ |V£| 2 /H 



(4.4) 



where C depends on 9 and a. Since J B u(.,t)£ A is bounded, and f. J B v p is bounded, we 
obtain an estimate of the gradient: 

(l + u) - 1 |Vu| 2 £ A S C. 




t J Bo 



Next recall the Gagliardo-Nirenberg estimate: let m > 1,7 G [l,+oo) and v G [0,1] such 
that 

1 1 1 , 1 - 1/ 
then there exists C = C(N, m, v, p)>0 such that for any w G W 1,2 (.Bp) n L m {Bp), 



\ w ~ M\li{b p ) <C\\\Vw\\\ v L 2 {Bp) \\w - HliUV) 



(4.5) 
(4.6) 



We apply it to w(x,t) = (1 + u(x,t))P , and 



a 1+a 2 2 

which satisfy (|4.5p . Therefore, for any t G (0,0) , 



m = -, 



(4.7) 



(l + n(.,t))^-^(t) 7 <C / (l + uM^lVuM)! 2 x 



Bp 



1//3 



(1-^)7/3 



Now IIW. 



L°°((0 6»)) — C because /3 G (0, 1) and J B u(-,t) is bounded; in turn we get 



Bo 



(l + u(x,t)f -w(t) 



VP 



dx < C, 
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Therefore, 

f (l + u(x,t))^ dx < C f (1 + u{.,t)) a ~ l \Vu(.,t)\ 2 dx + C. 

J Bp J Bp 

Integrating on (0, 9) we obtain 

pB p 

(1 + u(t))^ dx < C. 



Observing that j3^ = a + l + 2/N, and a is defined by (|4.3p we conclude to (|4.2p . ■ 
In order of proving Theorem 11.21 we show the following dichotomy property: 

Proposition 4.3 Assume p,q > 1. Let (it, v) be a positive solution of the system in £1 x 
(0, T) . Let xq G Q. Then the following alternative holds: 

(i) Either there exists a ball B(xq,p) C Q such that J Q T f B / XQ p ^{u q + v p ) < oo and two 
Radon measures and m2, p on B(xo,p), such that for any ip £ C®(B(xq, p)), 



lim / u(.,t)ip = / ipdmip, lim / v(.,t)tp = / ipdm,2 p , (4.8) 

t -*°JB(io,p) ^B(a:o,p) ' ^ J B(x ,p) J B(x ,p) 

(ii) Or for any ball B(xq, p) <zVt there holds §b(x p)( u9 vP ^ = 00 an ^ then 

lim/ (u(.,t) +v(.,t)) = oo. (4.9) 



Proof, (i) Assume that there exists a ball B(xq, /))cfi such that Jq Jb( Xo p)( u<1 + yP ) < °°- 
Then (j4.8|) follows from Lemma l4.1i 



(ii) Suppose that for any ball f Q f } 



B(x ,p) 



U q + 



oo. Consider a fixed p > such 



that B(xo,p). We can assume xq = 0. We choose the test function £ A , where £ is defined at 
(I3T3D and A > 2max(p',g / )- Then 



As above from ()3.7p . since A is large enough, 



i; p £ 



u(A£ ; 



/« 


A£ A 




I B p 




J B 



where C depends on p. Let < t < 9 < T. Consider X, Y, Z, W defined by (|4.ip . Then we 
find with new constants C > 

X t it) + w(t) s cz x i\t) s^ + c, 

Y t (t) + Z(t) S CW l l p (t) <, ^ + C. 
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By addition 

Z + W 

(X + Y) t (t) + ^-(t)SC 
By hypothesis Z + W $ L 1 ((0,T)), then 



lim(X(t)+y(t)) = lim / (u(.,t)+v(.,t))t = oc, 



B„ 



thus 



lim / (u(.,t) + «(.,t)) = oo (4.10) 



for any p < p, and the conclusion follows, since p is arbitrary. ■ 
As a direct consequence we deduce Theorem 11.21 

Proof of Theorem 11.21 Let 

(«(., t) + «(., t) < oo L, 



^ = j^o G ^ : 3p > 0, B(xo,p) C f2, lim sup ^ 



B(x ,p) 

and S = Q\H. Then 7£ is open, and from proposition 14.31 there exists unique Radon 
measures p±,p2 on 1Z such that (|1.7p holds, and (|4.9p implies (|1.8p on any open set U such 
that U n S / 0. ■ 

Next we give more information when p, q are subcritical. 

Proposition 4.4 Assume < p,q < 1 + 2/ N. Let (u,v) be a positive solution of the system 
in 0, x (0, T) . Let xo € Q. Then then the eventuality (ii) of Theorem \4-3\ is equivalent to: 

(in) for any ball B{xq,p) C SI there holds 

u q = oo and f f v p = oo. ( 4 -ll) 

o Jb(x , p ) Jo Jb(x ,p) 

Proof. It is clear that (iii) implies (ii). Suppose that (hi) does not hold, and reduce to 
xq = 0. Then there exists a ball B p such that for example 



T 



v p < oo. 



Then for any p < p f B _ u(.,t) is bounded as t — > 0, from Lemma (j4.1|) . Since q < 1 + 2/iV, 
we obtain 



u q dx S C, 

for any p' < p, from Lemma 14.21 Then (ii) does not hold. 



t JB p , 
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Remark 4.5 Under the assumption (ii) or (Hi) of Proposition \4-4\ f or an V ball B p = 
B(xq,p) C ft, j B u(.,t) and J B v(.,t) are unbounded near 0, from Lemma \4-S\ But we 

cannot prove that lim^o Jb(x p) u (->t) = oo or lim^o Jb(x p) v (-^) = °°> even in the case 
p, q > 1 where \4-ty holds. 



We give a last result concerning the case where the two equations are sublinear. 



Proposition 4.6 Assume < p, q 5= 1. Let (u,v) be a positive solution of the system in 
ft x (0,T) . Then there exist two nonnegative Radon measures /J>i,fJ<2 on ft, such that for 
any </> G C°(fi) 

lim / u(.,t)i/j = / ipdfix, nm / = / 4>d(i2, 

Proof. Consider any ball B(xo,p) C ft, and assume xq = 0. Consider again X,Y,Z,W 
defined by (14. ip . Here we find 

W(t)= [ ^(.,t)£ A ^/ (v(.,t) + l)ti x SY(t) + C, 

J Bp J Bp 

z(t)= [ u i(.,t)es[ (u(.,t)+i)esx(t)+c, 

J Bp J Bp 

and 

— {e XXl ^X{t)) + e XXl ^W(t) ^ 0, —(e^^Yit)) + e XXl ^Z(t) ^ 0, 
at dt 

then the function <£> = e AAl ^'(A"(t) + Y(t)) satisfies + <&(t) + Ce XXl ^ ^ 0, that is 

(e*($0) + C(l + AAi )P ) _1 e AAl 'P*) / ^ 0. Then $(t) has a limit as f — > 0. ■ 

Open problems: 

1) Can we extend Theorem 11.21 to the case pq > 1? 

2) Can we extend Proposition 14.61 to the case pq < 1? 



5 Removability results 

Here we prove the removability of punctual singularities when p and q are supercritical. 

Proof of Theorem 11.41 We can assume that q > p > 1 + 2/N. Let cj be a regular domain 
such that u CC fi\{0} and let T\ < T. Then from flTJJ]) u,v G L°°(0, Xi; L^u;)); then from 
LemmaEOl u G x (0,Ti)) and t> G L p (o; x (0,Ti)). As in Theorem 2], step 3, the 

functions defined on uj x (— T, T) by 

satisfy £t G Lf oc (w x (0,T)), € G Lf oc (w x (0,T)), and 

u t - Au + v p = 0, v t - Av + u q = 0, inP'(w x (-T,T)). (5.1) 
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It follows that u,v e C 2 ^(uj x [0,T)) and 

u(x,0) = v(x,0) = 0, Vi G oj. 
Since p < q, we have vP ^ vfl + 1 from the Young inequality, thus the function 

<? = 2 (1 - p)/p (n + u) 
satisfies g G L p (w x (0,Ti)), flr(ar; 0) = 0onw\ {0} and 

gt-Ag + cf^l 

in w x (0, T). Following [6j Theorem 2], step 4, we deduce the key point estimate: there 
exists C = C(N,p) and p > such that -E>(0, 2p) C $7, and T\ <T such that 

t} = (t+ | x |^i/ (p -i) + C ' V ( x ' *) G P) x (°> Tl ) • ( 5 - 2 ) 

Since p > 1 + 2/N, it implies that g € L 1 (i?(0, p) x (0,Ti)). From [61 Theorem 2], step 5, 
it follows that g £ L p (B(0, p) x (0,Ti)), thus also u and u. We claim that a better estimate 
holds, adapted to the system: 

1 [ v p <oo and / 1 / u q < oo. (5.3) 
o JB(o,p) Jo Jb(q,p) 

Indeed, consider a function (gP((1x (— T, T)) with values in [0, 1], such that £ = 1 on 
i?(0, p) x (0, Ti), and a function x £ C 00 ^), nondecreasing, with x{t) = for t ^ 1, x(i) = 1 
for i 2l 2; let Xfc(^) = x(^) f° r an y & > 1- Setting 

D fe = j( x ,i) : < \x\ 2 + t< 2/A:} , 

and using the test function 

Pkfat) = Xk{\x\ 2 + t)((x,t), 

we obtain 

/ u q <Pk = / v(ip k ) t + / / vA<p k , 

JBp JO JB P JO JB P 

/ l^fc = / / U(lfk)t + / / uA(f k . 

jb p Jo Jb p Jo Jb p 

Consequently 

f [ u q ip k SCk [ f v + C, [ [ v p ip k SCk [[ u + C. (5.4) 

Jo J B p JJo k Jo J B p J J D k 

Next from (|5.2|) . we have 

(B+B)g ^ WiW-> +c H (5 - 5) 



D k JJD 
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Hence (|5,3p follows from (|5.4p . (|5.5|) and the Fatou Lemma. As a consequence of (|5.3p . 
u e L\ oc (n x (-T, T)) and 5 € Lf oc (0 x (-T, T)) . 
Following [6, Theorem 2], step 6, we have 

(« + «) (Kxfe)tCI + l(Axfc)Cl + |Vx*| |VC|) ^ Cfc // (« + «), 
and from the Holder inequality 

fc/7" (u + v)^k( [f (u + vy) 1 Vfcl W ^cf// (« + «FV P - ( 5 - 6 ) 

iifl t \JJD k J \JJD k J 

Since the right hand side of (|5.6|) tends to from (|5.3|) . we can pass to the limit as — >• oo 
in (??), and obtain 

T [ u%= F f v( t + F f vA(, [ T [ v*>( = F ! u( t + [ T [ nAC, 
o Jb p Jo Jb p Jo Jb p Jo Jb p Jo Jb p Jo Jb p 

and then 

ut - Au + i p = 0, v t - Av + u q = 0, in V (ft x (-T, T)) . 

Therefore u,v G C 2 ' 1 (fi x (-T, T)), and n(x,0) = u(z,0) = on O. ■ 

Open problem: In the elliptic problem (jl.4p in i?(0, 1)\{0}, it was shown in [H 
Corollary 1.2] that the singularities at are removable as soon as 

max(2a, 26) ^ N - 2. 

In the case of system (ll.ip , an open question is to know if the initial punctual singularities 
at are removable whenever 

( h \< N 
max(a, o) S — , 

a condition which is obviously satisfied when p, q ^ 1 + 2/N. 
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